AD-A0J1  432  BATTELLE  COLUMBUS  LABS  OHIO  F/q  11/7 

THREE  OIMENSIONAL  STRESS  ANALYSIS  OF  A LAMINATED  PLATE  CONTAINI--ETC ( 
APR  76  E F RYBICKX'  0 « SCHMUESER  F33615-7«-C-5096 

UNCLASSIFIED AFML-TR-76-32  . NL 

• sl0VHI|HWnHP 


A 031432 


ADA031432 


THREE-DIMENSIONAL  STRESS  ANALYSIS  OF  A LAMINATED 
PLATE  CONTAINING  AN  ELLIPTICAL  CAVITY 


BATTELLE 

COLUMBUS  LABORATORIES 
505  KING  AVENUE 
COLUMBUS,  OHIO  43201 


APRIL  1976 


FINAL  REPORT  FOR  PERIOD  DECEMBER  1973  DECEMBER  1974 


Approved  for  public  release;  distribution  unlimited  j 

" t 


AIR  FORCE  MATERIALS  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT  PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


D D C 

. . nr 


NOV  1 ,9T6 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a definitely  related  Govern- 
ment procurement  operation,  the  United  States  Government  thereby 
incurs  no  responsibility  nor  any  obligation  whatsoever;  and  the  fact 
that  the  Government  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to 
be  regarded  by  implication  or  otherwise  as  in  any  manner  licensing 
the  holder  or  any  other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use,  or  sell  any  patented  invention 
that  may  be  related  thereto  in  any  way. 

This  report  has  been  reviewed  and  cleared  for  open  publication  and/or 
public  release  by  the  appropriate  Office  of  Information  (01)  In 
accordance  with  AFR  190-17  and  DODD  5230.9.  There  is  no  objection  to 
unlimited  distribution  of  this  report  to  the  public  at  large,  or  by 
DDC  to  the  National  Technical  Information  Service  (HTIS). 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


N^TPaga^o  

Project  Engineer 

FOR  THE  DIRECTOR 
S.  W.  Tsai,  Chief 

Mechanics  and  Surface  Interactions  Branch 
Nonmetalllc  Materials  Division 


Copies  of  this  report  should  not  be  returned  unless  return  Is  required 
by  security  considerations,  contractual  obligations,  or  notice  on  a 
specific  document. 


»'•  'Wl  - IT  VtPTfffhft  76-700 


UNCLASSIFIED 

security  classification  of  this  rage  (Whmn  p«i«  gnitnuj 

[ RETORT  documentation  page 


KTR-76-32 


12  GOVT  ACCESSION  NO. 


IT  title  (M 


/THREE  DIMENSIONAL  STRESS  ANALYSIS  OF  A / 

Sj  / LAMINATED  PLATE  CONTAINING  AN  ELLIPTICAL  j t/ 
— 1 / CAVITY * ' 


READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

3 recipient's  catalog  number 


/Finale-Report  • 

Dec  B73  tm-Dec  <974, 


J Edmund  F.  /iRyblckl^BF  David  W.  Schmueser 


F33615-7Hf5096 


* PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Bat telle,  Columbus  Laboratories 

505  King  Avenue  y, 

Columbus,  Ohio  43201 

II  CONTROLLING  OEEICE  NAME  AND  AOORESS 

Air  Force  Materials  Laboratory  (AFML/MBM)  y 

Air  Force  Wright  Aeronautical  Laboratories  — - 

Wright-Patterson  AFB,  Ohio  45433 

u MONITORING  AGENCY  name  S ADDRESS«f  dllfarani  Iron*  Controlling  Olflca) 


10  RROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 


Project  7340 
Task  734003 


' Apr£l  1976 

"U  NUMBER  Of  PAGES 

. -42.-..- u 

15  SECURITY  CLAS$>f 

Unclassified 


15m  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


Ms  Distribution  statement  (of  thia  Report) 


Approved  for  public  release;  distribution  unlimited. 


I *7.  DISTRIBUTION  STATEMENT  (of  tho  abatract  anlarad  In  Block  20,  It  dltfaranl  from  Roporl) 


»•.  supplementary  notes 


It.  KEY  WOROS  (Coollnue  on  reran*  alda  If  nocaaaary  and  Idanllty  by  block  number) 

Stress  Analysis  Delamlnatlon 

Elliptical  Hole  Laminated  Plates 

Circular  Hole 
Crack 

Interlaminar  Stresses 

20  ABSTRACT  (Conllnua  on  rararaa  alda  If  nocaaaary  and  Idanllty  by  block  number) 

— - An  analysis  procedure  is  developed  to  study  the  stress  distribution  In  a 
laminated  plate  containing  an  elliptical  cavity  under  tenalde  loading  and  a 
uniform  temperature  change.  Also  considered  is  the  stress  field  around  a 
radial  crack  emanating  from  a hole  in  a laminate.  The  analysis  follows  from 
an  extension  of  a previously  developed  finite  element  approach. 


DD 


EDITION  OP  I NOV  «»  IS  OBSOLETE 


UNCLASSIFIED  ^ fy 

SECURITY  CLASSIFICATION  OP  THIS  PAGE  I 


' ■****<>', -r. 


PREFACE 


In  this  work,  a previously  developed  solution,  based  on  a higher 
order  finite  element  simulation,  is  extended  to  treat  the  problem  of  a 
composite  laminate  containing  an  elliptical  cavity.  Loading  conditions 
consist  of  axial  tension  plus  a uniform  temperature  change.  This  class 
of  problems  is  fundamental  to  an  understanding  of  local  failure  mechanisms 
which  prevail  in  the  vicinity  of  cutouts  or  cracks  in  composite  laminates, 
where  the  influence  of  lamination  geometry  must  be  recognized.  Furthermore, 
the  inclusion  of  thermal  stresses  in  the  analysis  allows  the  calculation 
of  curing  stresses  induced  in  the  fabrication  process.  Computer  program 
listing  and  input  data  instructions  for  this  stress  analysis  routine  may 
be  obtained  from  Dr.  N.  J.  Pagano  (AFML/MBM) , the  project  monitor  for 
this  program. 
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SECTION  I 
INTRODUCTION 


One  of  the  unique  characteristics  of  laminated  composite  materials 

is  the  delamination  mode  of  failure  that  occurs  at  free  edges.  This  behavior 

has  been  investigated  experimentally  and  by  mathematical  stress  analysis 

models  for  laminates  with  straight  free  edges.  Experimental  studies  by 

Pipes,  Kaminski,  and  Pagano1  have  shown  that  the  laminate  stacking 

sequence  can  affect  the  static  strength  of  laminates.  Similar  effects  of 

[2] 

stacking  sequence  were  found  by  Foye  and  Baker  for  fatigue  loading 

conditions.  Mathematical  stress  analyses  of  straight  free  edges  of  laminate 

plates  have  shown  that  changes  in  stacking  sequence  can  change  the  behavior 

[ 3 4] 

of  the  interlaminar  stresses.  In  particular,  investigators  * have  found 

that  the  sign  of  the  transverse  normal  stress  oz  changes  from  tension  to 

compression  as  the  stacking  sequence  is  changed.  Furthermore,  Pagano  and 

Pipes^’^  have  shown  that  high  tensile  stresses  are  associated  with  the 

decreased  laminate  strengths  reported  in  References  [1]  and  [2].  These 

observations  point  to  the  importance  of  understanding  the  stress  behavior 

of  laminates  near  a free  edge.  Through  such  an  understanding,  Lackman  and 

[ 8] 

Pagano  have  successfully  reduced  interlaminar  stresses  at  the  free  edge 

to  alleviate  the  effects  of  the  delamination  mode  of  behavior. 

While  the  straight  free  edge  can  be  studied  as  a problem  with 

two-dimensional  variations  in  stresses  and  displacements,  the  curved  free 

edge  is  inherently  a three-dimensional  problem.  Thus,  the  circular  hole 

problem  is  more  difficult  to  treat  in  terms  of  a stress  analysis  than  the 

f 91 

straight  free  edge  problem.  Dana  and  Barker  have  investigated  the 

three-dimensional  stress  analysis  of  four-ply  (0/90)  and  (45/-45)s 

laminates  with  circular  holes  and  found  that  the  sign  and  magnitude  of  the 

oz  distribution  around  the  hole  can  change  with  changes  in  stacking  sequence. 

Rybicki  and  Hopper^  ^ have  developed  an  analysis  for  the  interlaminar  stress 

distributions  for  a six-ply  symmetric  laminate  containing  a circular  hole. 

Experimental  investigations  by  Daniel,  Rowlands  and  Whiteside^^  on  the 

effects  of  stacking  sequence  on  the  strength  of  laminated  plates  with 

holes  Indicated  that  variations  in  stacking  sequence  effected  the  laminate 

f 121 

strength  and  failure  mode.  Other  studies  for  example  by  Greszczuk  , 
Waddoups,  et  al.^^,  and  Waszczuk  and  Cruse^^  have  predicted  failure 
of  a composite  plate  containing  a hole  and  made  comparisons  with  experimental 
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data.  These  studies  are  a first  step  toward  understanding  the  behavior  of 
laminates  containing  a circular  hole.  The  stress  analyses  for  these  studies 
did  not  include  interlaminar  stresses.  Two  basic  three-dimensional  problems 
that  can  provide  insight  into  understanding  laminate  behavior  near  a curved 
free  edge  are  considered  here.  One  problem  is  the  state  of  stress  that 
exists  around  a circular  hole.  The  other  problem  is  the  behavior  around 
a noncircular  hole  such  as  an  elliptical  hole. 

In  the  following,  a description  of  an  analysis  procedure  to 
treat  stresses  around  a circular  or  elliptical  hole  in  a laminated  plate 
due  to  a uniform  loading  distant  to  the  hole  is  described.  The  analysis 
procedures  have  been  implemented  into  an  existing  three-dimensional  finite 
element  stress  analysis^ ^ and  numerical  results  for  example  cases  are 
presented. 

SECTION  II 

METHOD  OF  STRESS  ANALYSIS 

The  method  of  stress  analysis  is  based  on  the  equilibrium  finite 
element  analysis  described  in  Reference  [10],  The  details  of  the  analysis 
are  given  in  that  reference.  A summary  of  the  analysis  is  given  here  as 
background  information  and  extensions  are  described  in  detail. 

In  the  analysis,  each  ply  is  modeled  as  a homogeneous  orthotropic 
material  with  linear  stress-strain  behavior.  Symmetric  laminates  are 
considered  with  applied  stress  loading  conditions  distant  to  the  hole. 

The  analysis  is  based  on  a three-dimensional  equilibrium  finite  element 
representation  for  the  stresses.  Three  Maxwell  stress  functions  are  associated 
with  each  element  to  satisfy  equilibrium  exactly  within  each  element. 
Equilibrium  between  contiguous  elements  is  satisfied  exactly  by  equating 
unknown  coefficients  between  elements.  All  stress  boundary  conditions  are 
satisfied  exactly.  In  the  equilibrium  approach,  displacement  boundary 
conditions  and  compatibility  are  satisfied  approximately. 

The  following  sections  contain  a description  of  the  extensions 
of  the  existing  analysis  to  treat  an  elliptical  hole  in  a laminated  plate. 

The  results  of  interest  for  the  elliptical  hole  problem  are  the  interlaminar 
stresses. 
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In  Reference  [15],  Ryblcki  and  Hopper  examined  an  approach  to 
transform  the  problem  of  a plate  containing  an  elliptical  hole  into  a 
problem  of  a plate  containing  a circular  hole.  This  is  accomplished  by 
starting  with  the  conformal  transformation  for  mapping  an  ellipse  into  a 
circle.  The  transformation  is  applied  to  the  governing  equations  of  the 
elliptical  hole  problem.  The  result  is  a description  of  an  equivalent 
circular  hole  problem  with  different  material  properties  and  loading 
conditions.  A transformation  of  the  circular  hole  solution,  back  to 
the  elliptical  hole  problem,  is  obtained  based  on  the  inverse  of  the 
original  transformation.  This  approach  was  extended  to  the  three-dimensional 
case.  The  procedure  is  given  in  the  following. 

A typical  lamina  is  shown  in  Figure  1,  The  stress-strain  relations 
for  each  lamina  of  the  laminate  have  the  form 
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where  the  S^j's  are  the  compliance  coefficients,  a's  and  T's  are  the  stresses, 
and  or^AT's  are  the  strains  due  to  temperature  change.  The  stresses  are 
related  to  the  Maxwell  stress  functions  ^2'  an<*  *3 
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The  first  step  is  to  express  the  strains  of  Equation  (1)  in  terms 
of  the  stress  functions,  xj,  by  combining  Equations  (2)  and  (1).  Next,  the 
compatibility  equations  are  expressed  in  terms  of  the  stress  functions. 

The  transformation 


was  applied  to  all  the  equations.  One  axis  of  the  ellipse  lies  along  the 
x-axis.  This  axis  has  a length  of  1.0.  The  value  of  £ is  then  the  length 
of  the  axis  of  the  ellipse  in  the  y-direction.  The  (x'  ,y'  ,z‘ ) system 
pertains  to  the  equivalent  laminate  with  a circular  hole  and  different 
material  properties.  For  brevity,  only  the  results  of  the  transformation 
are  listed.  The  material  properties  for  the  circular  hole  problem  (S^j 
and  o^)  are  related  to  those  for  the  elliptical  hole  problem  (S^  and 
and  to  0 as  shown  in  the  following  equations. 

S'  - Su/B4 


s;2  - S12/e 

S13  " S13/B 
S16  “ S16/P 
S23  “ S23/0 
S26  ' S26/B 
*1  “ S33  ■ S33^P 

<*l  e a2  S36  " S36/03 

»3  * »3/02 

»«  - V9  S45  - S45/03 
S55  * S53/02 
SG6  " S6G/P 

The  stress  functions  are  transformed  as  follows. 
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(5) 


X3  = X3 

where  the  primes  denote  the  circular  hole  problem. 

The  relationship  between  the  stresses  for  the  circular  hole 
problem  and  the  elliptical  hole  problem  are 
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where  the  primed  stresses  refer  to  the  circular  hole  problem  and  the  unprimed  t 
stresses  refer  to  the  elliptical  hole  problem. 

The  stress  boundary  conditions  around  the  elliptical  hole  were  j 

transformed  to  boundary  conditions  for  the  equivalent  circular  hole  problem. 

The  algebra  for  this  step  is  given  in  Appendix  A while  Appendix  B summarizes  | 
the  constraint  equations  to  satisfy  the  stress-free  conditions  at  the 
circular  hole. 

The  procedure  for  including  these  transformations  into  the 
computer  program  starts  with  the  elliptical  geometry,  the  material  properties 
and  the  loading  conditions.  Through  Equations  (4)  and  (6)  the  program 
transforms  the  properties  and  stress  boundary  conditions  to  those  of  an 
equivalent  circular  hole  problem.  After  the  solution  is  obtained,  the 
stresses  for  the  equivalent  circular  hole  problem  are  evaluated  in  the 
rectangular  (x' ,y‘ ,z‘)  coordinate  system.  These  stresses  are  transformed 
to  the  stresses  for  the  elliptical  hole  in  the  (x,y,z)  coordinate  system  by 
using  Equations  (6).  Additional  transformations  of  the  Mohr’s  circle  type 
were  included  in  the  program  to  obtain  the  tangential  and  interlaminar 
stresses  around  the  elliptical  hole. 
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The  stress  analysis  has  the  flexibility  to  include  thermal  stresses 
and  also  to  evaluate  energy  release  rates  J>.  This  can  be  done  by  solving 
the  problems  with  different  flaw  lengths  and  evaluating  from  the  change 
in  stored  energy  of  the  system  per  unit  area  of  extension  in  flaw.  The 
procedures  for  doing  the  thermal  stress  analysis  and  the  energy  release 
rate  calculations  are  outlined  in  Appendices  C,  D,  and  E.  Appendix  F contains 
description  of  the  input  data  for  the  computer  program. 

SECTION  III 
RESULTS 


Several  problems  were  run  wihh  the  equilibrium  finite  element  stress 
analysis  program.  These  problems  involve  a four-ply  (90/0) g laminate.  The 
configurations  and  loading  conditions  were  a laminate  containing  a circular 
hole  under  inplane  loading  and  a laminate  containing  an  elliptical  hole 
also  with  inplane  loading. 

The  principal  material  properties  for  the  plys  are 

Ex  = 30.  E2  = E3  = 3.0 


The  radius  of  the  hole  was  1.0.  The  thickness  of  each  ply  was 
0.2  inch.  Stress  boundary  and  conditions  were  applied  at  a radius  of 
12  inches  from  the  center  of  the  hole.  These  stresses  were  evaluated 
from  a laminated  plate  theory  representation  for  a uniform  resultant  stress 
in  one  direction. 

The  first  configuration  was  a (90/0)  laminate  containing  a circular 
hole  and  lnplane  loading  in  the  x-direction.  A reference  solution  for  this 
problem  was  obtained  by  three-dimensional  displacement  finite-element 
computer  program  called  SAP  IV.  The  obtained  solution  and  the  reference 
solution  for  the  interlaminar  stress,  oz,  around  the  circular  free  edge 
are  shown  in  Figure  2.  Figure  3 shows  the  distributions  along  the 
x and  y directions  for  the  equilibrium  finite-element  analysis. 


FIGURE  2.  DISTRIBUTION  OF  INTERLAMINAR  STRESS,  <rt  , 

AROUND  A CIRCULAR  HOLE  IN  A (90/0),  LAMINATE 
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Next  the  problem  of  an  elliptical  hole  in  (90/0^  laminate  was 
considered.  The  ratio  of  the  minor  axis  to  the  major  axis  was  0.5.  The 
loading  condition  was  inplane  loading  parallel  to  the  minor  axis  as  shown 
in  Figure  4.  The  interlaminar  stress  distribution  for  oz  at  the  midsurface 
is  also  shown  in  Figure  4. 

SECTION  IV 

SUMMARY  AND  DISCUSSION 

The  problem  of  the  three-dimensional  stress  behavior  around  holes 
in  laminated  plates  is  considered.  This  problem  is  inherently  three- 
dimensional  and,  as  discussed  in  Reference  [10],  the  characteristics  of  the 
interlaminar  stress,  o2,  cannot  be  inferred  from  equilibrium  considerations 
alone  as  has  been  shown  to  be  possible  for  the  straight  free  edge  problem. 
Experimental  studies  and  mathematical  stress  analysis  studies  addressing  the 
free  edge  problem  are  found  in  the  literature.  Experimental  studies  show 
that  stacking  sequence  has  a pronounced  effect  on  laminate  strength  and 
fatigue  life,  Indicating  a dependence  on  the  sign  of  the  interlaminar  stresses. 
However,  there  is  a need  for  a better  understanding  of  the  basic  behavior 
of  the  interlaminar  stress  distributions  around  holes  in  laminated  plates. 

Two  problems  that  can  provide  Insight  into  the  behavior  of  the 
Interlaminar  stress  distribution  near  a curved  free  edge  are  considered 
here.  One  problem  is  to  determine  the  stresses  around  a circular  hole.  The 
other  problem  pertains  to  a laminate  containing  an  elliptical  hole.  A mathe- 
matical stress  analysis  procedure,  based  on  the  three-dimensional  finite 
element  equilibrium  model  of  Reference  [10],  was  developed  to  handle  the 
problems  of  an  elliptical  hole. 

The  first  problem  was  that  of  a circular  hole  in  a (90/0)  laminate 
loaded  with  a uniform  inplane  extension  in  the  direction  of  the  x-axls. 

Here  a reference  solution  was  available.  The  reference  solution  was  obtained 
by  a three-dimensional  compatible  finite  element  stress  analysis  called 
SAP  IV.  The  comparison  of  the  interlaminar  stress  distributions  around 
the  circular  hole  was  good.  The  equilibrium  element  stress  analysis  gave 
higher  stress  gradients  and  higher  peak  values  of  oz,  but  both  analyses 
revealed  a changing  sign  of  oz  around  the  circular  edge. 
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The  next  problem  was  that  of  (90/0)g  laminated  plate  containing 

an  elliptical  hole  with  a ratic  of  minor  axis  to  major  axis  of  1/2.  The 

loading  was  in  the  direction  parallel  to  the  minor  axis  or  in  the  y direction. 

No  reference  solution  was  available,  but  the  stress  distribution  for  a 

\ 

around  the  elliptical  hole  appears  reasonable  if  one  considers  this  case  is 
similar  to  loading  a (0/90) s laminate  containing  a circular  hole  with 
inplane  loading  in  the  x direction. 

The  agreement  between  the  solution  generated  here  and  the  reference 
solution  gives  a degree  of  confidence  that  the  computet  program  for  the 
analysis  is  in  working  order.  This  program  can  be  used  to  generate  a 
catalogue  of  solutions  for  various  lay-up  angles  and  stacking  sequences 
to  provide  a basis  for  a combined  analyt ical /experimental  program  to  obtain 
a better  understanding  of  the  interlaminar  stress  and  strain  behavior  in 
laminates  containing  holes. 


APPENDIX  A 


TRANSFORMATION  OF  BOUNDARY  CONDITIONS  BETWEEN 
THE  CIRCULAR  AND  ELLIPTICAL  HOLE  BOUNDARY 


The  transformation  of  stress  boundary  conditions  between  the  circular 
and  elliptical  hole  problems  can  be  represented  by 


faE) 


Nez 


[c]  [ac] 


NQZ 


(A- 1) 


where 


{aEl 


NQZ 


E 'l 

C ' 

aN 

aN 

E 

C 

CTQ 

ctq 

E 

c 

az 

r c. 

CTz 

E 

V 5 ^NQZ  * < 

c 

tnq 

TNQ 

E 

c 

tnz 

tnz 

E 

c 

. T0Z 

. T0Z  j 

and  Ic]  is  a transformation  matrix. 

Let  [ T ] represent  the  transformation  matrix  which  relates  stresses 
in  polar  coordinates  to  stresses  in  cartesian  coordinates.  For  the  circular 
hole , 


r ci 

{<*  1 


NQZ 


{TC1  f<TC] 


XYZ 


(A -2) 


and  for  the  elliptical  hole, 


f Ei 

1°  1 


NQZ 


[TC1 


XYZ 


(A-3) 


The  matrix  ( C I in  Equation  ( A— 1 ) can  be  derived  from  Equation  (A-l), 
Equation  A-3,  and  the  transformations  which  relate  cartesian  stresses  of 
the  circular  hole  problem  to  cartesian  stresses  of  the  elliptical  hole 
problem.  These  trans format  ions  are 
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(A-4) 


E C . 2 
CTX  * VB 


Y 

E 


C , 2 
CT  Z/8 


XY 

C 

rxz 

E 

rYZ 


TXY/B 


c /Q  2 
TXZ/B 


YZ 


where  0 is  defined  by  the  cartesian  coordinate  transformation, 


X 


0Y' 

„C 


(A -5) 


From  Equation  (A-2), 


faCl 


XYZ 


[tC]-1£«jC3 


Y0Z 


(A-6) 


where 


„C  ,„2 

T /B 


[I0']-1 


tc3/b2 

T?/8 

Vs 


(4-7) 


C t h C ™ 1 

and  denotes  the  i row  of  (T  ] 

Using  Equation  (A-6)  and  Equation  (A-3), 


Thus, 


N0Z 


Ite]  [oEl 


XYZ 


[TE](TC  ] ~ 1 f CTC  } 


N0Z 


(A-8) 


(C) 


Ite](tc/]'1 


From  Equation  (A-8)  it  can  be  shown  that  the  boundary  conditions 

E E — E 

for  the  elliptical  hole  problem  are  oNOr,0,Z)  * TNZ(r,0,Z)  * TNgCr,0,Z)  •>  0 
when  the  boundary  conditions  for  the  circular  hole  problem  are  o^(7,0,Z)  * 

TN0(7'e,Z)  * tNZ(7»9*Z)  * 0 
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APPENDIX  B 


SATISFYING  STRESS  BOUNDARY  CONDITIONS  BY 
PRESCRIBING  APPROPRIATE  VALUES 

The  three-dimensional  finite  element  model  used  in  the  analysis 
of  the  laminated  plate  is  shown  in  Figure  B-l.  The  stress  boundary  conditions 
applied  to  this  model  are  listed  below. 


Table  B-l 


Locat ion 

Boundary  Conditions 

r " ro 

°r(x0,9,z)  m Tr0(ro»9‘Z)  * TrZ(rO»0,Z)  * ° 

r * rs 

Uniform  Stress  Field 

Z » 0 

T0Z(r,0,zo>  * T0r(r,e,ZO)  * 0 

z - Zs 

‘’zCr.e.Zg)  - T0Z(r,6,Zs)  * T0r(r,9,Zs)  « 0 

The  following  sections  describe  how  these  boundary  conditions  are  met  by 
prescribing  appropriate  F^*^  coefficients  which  appear  in  the  expressions 
for  the  stress  functions,  These  stress  functions  have  the  form 

X,  - its^V'VV1 

I j K 

where  I • l,  2,  or  3. 


Boundary  Conditions  at  R * R 


fl 


On  the  inner  hole  surface,  P^  * P^  * 1,  while  all  other  P^'s 
are  zero.  The  expressions  for  ar(r0),  Tr0(ro),  and  TrZ(r0)  and  therefore 

W • = * {FuiQ]Rk(1-co*2si>  + Fui,^to’29  + 

J * (B-l) 

+ FS^(Q;'/r2)Rk  + - 0 
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T_«(rf>)  * E E ' F^QjR^cos0  • sine  - F ^'QtR,' 'cos08  in0  + 


r9vl0 


j * 

<3> p(3) , 


,<2)0  r"c 

ijkQjRkC 


- + Flji<Qj/t2>RJ  - 0 


(B-2) 


TrZ(r0>  • =;  l • SI"»  * 


j k 


1 jk  j 


+ Fi jk(Qj/r)RkCOS0  ’ Sin0  " F3jkQjRkCOs20/  “ ° 

(3) 

Setting  the  terms  containing  F in  Equations  (B-l)  and  (B-2) 

i jk 


(B-3) 


to  zero  gives 


z I F;j(Q)*/.!).k  * I z iSjBjHi,  . o , 

i fSm3/c)\  + F!]k«;/r2)Rk  - ° 


( B-4) 


(B-5) 


Two  sets  of  constraint  equations  are  generated  by  using  symmetric 
and  antisymmetric  terms,  respectively.  For  example,  the  expansion  using 
symmetric  terms  y » 1,3,5,. ..,11  is 

-»4- 

* {-F5»+rm+F3s£}  c°*69  + {-Fm+F»i+F32i}  cos8e  + (B-6> 

* {■F39k+F3!il,k  + F39k}  ~l#.  * ' ° 

k * l, 2, 3,4 

Therefore,  the  set  of  constraint  equations  is 
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Fa>  + F(2) 

33k  31k 


_pd>  + F(l)  + f(2)  „ 0 

* 33k  35k  33k 

„pd)  + Fd)  + p<2>  , o 

r 35k  37k  35k 


Fd)  + pd)  + p(2)  „ o 

•F37k  + F39k  + F37k  ° 


.pd)  + p(i)  + pd)  q 

39k  3,11, k 39k 

-F(1)  + F(2)  - 0 

3,ll,k  3,ll,k  * 


k - 1,2, 3,4 


The  expansion  using  antisymmetric  terms  y * 2, 4, 6,... 10  gives  the 
owing  set  of  constraint  equations. 

.pd)  + pd)  + p(2)  „ o 

* 34k  36k  34k  » 

pd)  + pd)  + p(2)  „ 0 

_F36k  + F38k  + F36k  ' ° * 


,Fd)  + pd)  + p(2)  „ o 

38k  3, 10,k  3,18,k  ’ 

pd)  + p(2)  , 0 

F3, 10, k F3, 10, k ° » 


k - 1,2, 3,4 


Boundary  Conditions  at  R « Rt 


The  following  conditions  satisfy  the  uniform  stress  boundary 


conditions  at  R « R[, 


a « <j  a m a and  t ■ t 

x x * y y xy  xy 


where  ?x,  oy,  and  Txy  are  prescribed.  By  letting 

1-2  2 1-2  2 - 2 
X3  * J °xRpsin  0+2  CTyRpcos  9 _ TxyFFco80slne  , 


the  uniform  stress  requirement  at  R ■ Rp  is  met  by  prescribing  the  following 

p(3)  , 

r ' values . 
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(3) 

« F(3> 

l - „2 
» - o R 

F/(3> 

- fill 

211 

212 

2 X F 

411 

412 

(3) 

231 

« f(3) 
232 

“ I Rj>y-x> 

f(3) 

431 

.f(3) 

432 

(3) 

. f(3) 

m R2 

fP> 

K F(3) 

221 

222 

xy  F 

421 

422 

a R 
x F 


RF(ay<rx) 


-2  TXy  Rp 


(B-9) 


f(3)  - f(3) 
2jk  bijk 


<■  0 


for  k m 3,4 


The  prescribed  values  of  F^,?  and  F^?  are  obtained  by  the  same 

ijk  ijk 

analysis  that  was  applied  at  the  inner  radius,  except  that  at  R • Rt,  ^2  " Pu  " ^ • 
while  all  other  P^'s  and  P f ' s are  zero.  At  R * Rt  the  analysis  gives  the 
following  prescribed  F^^  values  and  constraint  equations. 


F(1)  -F(2)  « 0 
ijk  f 2 jk 


F^  + p^2)  m 0 

*43k  41k  * 


_Ka>  + p(l)  + p(2)  „ 0 

*43k  45k  f43k  U • 

.F(l)  + Fd)  + F(2>  « o 

*45k  47k  45k  * 

,F(1)  + FU>  + p(2)  . o 

47k  49k  47k  * 

-F^  + F^  + F^  “ 0 
49k  4,ll,k  49k  * 

-F(1)  + F(2) 

r4,ll,k  4, ll,k 


for  all  j and  k , 


PF^1)  + f^2^  >=  0 

44k  42k  » 

_p^^  + p(D  + F^2^  * 0 

^ 44k  46k  44k 

_F(1>  + pd)  + p(2)  . o 

* 46k  48k  b46k  U * 

-F^  +F^  + F^2^  •=  0 

*48k  4,10,k  48k 

.p^1-)  + p(2)  „ 0 

4,10,k  4,10,k  » 

, k -=  1,2, 3, 4 


(B-10) 


These  equations  and  prescribed  F ! ^ values  apply  to  elements  2,4,  and  6 of 
the  finite  element  model. 


18 


Vi>  w. 


v 
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Boundary  Conditions  at  7.  * 0 

On  the  lower  plate  surface,  •=  1 while  all  other  R^'s  and 

R^'s  are  zero.  These  conditions  and  the  relationships 


p'  m -p' 
1 2 


p;-  -p* 


Q{  - Q["-  0 


(B-ll) 


give  the  following  prescribed  conditions  for  elements  1 and  2. 


r<D 

- F(l> 

F(1> 

113 

213 

ij3 

f(2) 

-<2) 

« r 

(2) 

F 1 1 3 

213 

i j 3 

0 

0 


j > 1 , i - 3,4, 5,6 


( B - 1 2 ) 


Boundary  Conditions  at  Z «=  Zs 

On  the  upper  plate  surface,  R2  *=  R,f  c 1 while  all  other  R^'s  and 
R^'s  are  zero.  These  conditir.s  and  Equation  (B-ll)  give  the  following 
prescribed  conditions  for  elements  5 and  6. 


(1) 

112 

"F(1)  , 

212  * 

'S  • ° ■ 

j 

1 , 

i 

- 3,4,5,6 

(D 

114 

m F^ 

214  * 

'15-  • 

j > 

1 , 

i 

* 3,4, 5, 6 

(2) 

112 

P(2) 

’ F212  • 

F!)2  ‘ 0 ■ 

j > 

1 , 

i 

K 3,4,5,6 

(2) 

114 

m F^2^ 

214  * 

'IS-  • 

J > 

1 . 

i 

* 3,4, 5,6 

( B-13) 


A comparison  of  the  appropriate  equations  developed  in  this  appendix 

shows  the  specified  F*.  values  to  be  compatible  in  elements  which  are 

1 JK 

bounded  by  more  than  one  surface. 
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APPENDIX  C 


THERMAL  STRESS  ANALYSIS 


This  appendix  describes  how  thermal  stresses  due  to  a constant 

uniform  temperature  were  included  in  the  laminated  plate  analysis.  The 

procedure  is  listed  in  terms  of  the  following  three  steps. 

Step  1.  Write  the  expression  for  the  part  of  the 

T 

complementary  energy  density,  uc , due  to 

thermal  strains  and  obtain  an  expression 

for  the  thermal  strains  in  the  (r,9,z) 

reference  coordinate  system. 

Step  2.  Express  u£  in  terms  of  the  stress  functions 

and  the  thermal  strains  of  Step  1. 

Step  3.  List  the  integrals  of  r,  9,  and  z that  are 

required  and  determine  how  to  evaluate  them. 

In  Step  1,  the  expression  for  the  complementary  energy  density 

T 

due  to  the  thermal  strains  is  denoted  by  uc  and  can  be  written  as 


T 

e a + 
r r 


T 

€ 

<TA 


T T T 

e a + y t + y 
z z 'rz  rz  Tqz 


T + V T 
qz  'rq  rq 


(C-l) 


The  material  is  assumed  to  have  three  principal  coefficients  of  thermal 
expansion  that  coincide  with  the  principal  material  directions.  These  are 
denoted  by  cr^,  , and  o^*  Figure  1 shows  the  principal  material  directions, 
the  (x,y,z)  and  (r,9,z)  coordinate  systems  for  a typical  lamina.  Expressions 
for  the  thermal  strains  in  the  (r,9,z)  systems  are  needed  for  the  analysis. 

To  obtain  these  expressions,  the  thermal  strains  are  first  transformed  from 
the  (1,2,3)  system  to  the  (x,y,z)  system  and  then  to  the  (r,9,z)  system. 

Transforming  the  thermal  strains  from  the  (1,2,3)  system  to  the 
(x,y,z)  system  gives 


T 

e 

x 

2 2 
■ (tt^cos  0 + c^siri  0) 

AT  , 

(C-2) 

T 

e 

y 

2 2 
■ (ojSin  0 + &2CC,S  0) 

AT  , 

(C-3) 
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T 


and 


ez  ■ q,3AT  • 

T 

VXy  * (2<y^cos  0 sin  0 - 2q'2cos  0 s in  0 ) AT  , 

T T . 
v “ V * 0 , 

’xz  ’ yz  ’ 


(C-4) 

(C-5) 

(C-6) 


where  0 is  shown  in  Figure  1 and  AT  is  the  temperature  change.  Since  0, 

At,  a^,  c^,  and  o/^  are  constants  for  any  lamina,  the  strains  in  Equations  (C-2) 
through  (C-6)  are  also  constants  for  any  lamina. 

The  result  of  transforming  Equations  (C-2)  through  (C-6)  to  the 
(r,9,z)  directions  is 


T 

T 2 

T 2 t 

(C-7) 

er 

» e m 

X 

•V 

e n + v mn 
y Txy  * 

T 

T 2 

T 2 T 

% 

* e n 

X 

+ 

e m - y mn 
y Txy  * 

(C-8) 

ez  " ez  » <C-9) 

T T T 2 2 T 

y * - 2e  mn  + 2e  mn  + (m  -n  ) v , (C-10) 

Yrq  x y Txy 


and 

T T 

V • v.  * o 
”rz  'Qz  * 

(C-ll) 

where 

m « cos  0 

(C-12) 

and 

n - sin  0 

(C- 13) 

The  strains  in  Equations  (C-7)  through  (C-ll)  are  functions  of  9.  Combining 
Equations  (C-7)  through  (C-ll)  with  Equation  (C-l)  gives  an  expression  for 
che  complementary  energy  due  to  the  thermal  strains  in  terms  of  stresses 
and  strains  in  the  polar  (n.9,z)  coordinate  system. 

In  Step  2,  the  stresses  that  are  needed  to  evaluate  uj  in  Equation 
(C-l)  are  or,  Oq,  and  The  stresses  T and  are  not  needed  because 

yjz  = * 0 from  Equation  (C-ll).  The  expression  for  u£  in  terms  of  the 

stress  functions  and  the  thermal  strains  is  obtained  by  expressing  the  stresses 
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in  Equation  (C-l)  in  terms  of  the  stress  functions  and  substituting  Equations 
(C-7)  through  (C-13)  into  Equation  (C-l).  The  result  for  a lamina  is 

T T 2 T 2 T 

u « (e  cos  0 ♦-  e sin  Q + y cos  9 sin  9)  x 
c x y xy 

ZEE  j*  P . (r)  Q.(o)  R"(z)  + 
ijk  1 ljk  1 J k 


- pi(r)  Qj(9)  Rk(z)  cos20  + 


+ P. (r)  Q (9)  R^(z)  cos2e  + 


+ E(3;  (P,(r)/r2)Q"(0)  R. (z)  + 


ijk  ' i 


+ F<3>  (P;(r)/r)  Q.(0)  Rk(z)}  + 

/ T 2 T 2 T 

+ 1 e sin  9 + e cos  a - v cos  a sin  9)  x 
\ x y xy 


??  Hik  Pi(r)  «j<9>  Rk<*>  'O*'9  + 

+ FS5k  pi(r)  Qj(e)  Rk(z)  + 


• FS  pi(r)  Q)(9)  RkW  'os«  + 


* Fuk  pi<t>  V9)  V*>}  + ‘I  - 


^ Wik  riM  ve)  v*>  + 

‘ Fi]k  pi(r)  Qj(0)  Rk(z)  cos2e  + 
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+ 2F[j^  (P|(r)/r)  Q'.  (0)  RR(z)  cos  0 sin©  + 


+ F<1>  (Pi(r)/r2)  Q'j’(e)  Rr(z)  cos2©  + 


+ (P|(r)/r)  Q^(ft)  RR(z)  cos2©  + 

- 2F[jk  (Pi(r)/r2)  Q!(©)  Rr(z)  cos  ©sin  © + 
+ F*2£  P'{(r)  Qj  (P)  Rk(z)  cos2©  + 


- 2F^2r  (P|(r)/r)  Qj(q)  RR(z)  cos  0sinfl+ 


+ (Pi(r)/r2)  Qj (0)  Rr(z)  + 


Fijk  (Pi<r)/r2)  Q'jO)  Rk(z)  cos2©  + 


+ (P^(r)/r)  Q.(0)  Rr(z)  + 


Fijk  Qj(9>  Rk(z)  cos2©  + 


(2)  2 

+ 2F  J jk  (Pt(r)/r  ) Qj(ft)  RR(z)  cos  © sin qj  + 


/ T T 2 T 

+ ^-2ex  cos  o sin  0+  2e^  cos  © sin  a - (2  cos  9_l)yxyJ  x 


r fl) 

UX  ^FjJk  Pjr)  QjCo)  RR(z)  ccs  © sin  © + 
i Jk 

(2) 

Ftjk  Pi(r)  Qj(<’)  Rk(z)  cos  ft  sin  © + 


Fijk  (Pi(r)/r>  Q](6)  Rk(z)  + 


+ FJjk  (Pt(r)/r2)  QJ(Q)  Rk(z)j 
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v i #; 


(C-1A) 


1 


r 


In  Equation  (C-14),  P^(r),  Qj(6),  and  Rj^(z)  are  the  assumed  functions  for 

the  finite  element  representation  of  the  three  stress  functions.  The 

F^‘?'s  are  the  unknown  coefficients  for  the  stress  functions  associated 
ijk 

with  P (r),  Q (a),  R.(z)  and  the  Lth  stress  function.  In  solving  thermal 

J a rrrr  t,  ^ 


stress  problems,  the  value  of 


dF 


(L) 

ijk  V 


u^dv  | must  be  evaluated  and 


subtracted  from  the  right-hand  side  of  the  system  of  equations.  It  can 
be  seen  from  Equation  (C-14)  that  the  following  integrals  are  needed  to 
evaluate  these  terms. 


| pi(r)rdr  , 


*1 

r2 


P (r)  f P (r) 

— 5 — rdr  * J dr 


‘ Pj(r) 


— rdr  - Pi(r2)  - , 


P'^(r)rdr  - r2P^(r2>  - - P^r^  + P.^)  , 

1 


J Qj(fl)f(e)db  , 


I 2 4 

where  f(0)  ■ 1 cos  0,  sin  0,  cos  0,  cos  0 sin  0, 

3 

cos  0 sin  0 , cos  0 sin  0 , 

42 

[ Q](e)g(e)do  , 

9, 


7i 


2 2 

where  g(0)  ■ 1,  cos  0 sin  0,  cos  0,  sin  0, 
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■ v t 


i‘q''(q)h(o)do  , 

J J 

91 

2 2 

where  h(q)  ■ 1,  cos  q sin  o , cos  q,  sin  0 , 

z2 

J Rk(z)dz  , 

Zl 


The  integrals  in  r can  be  evaluated  by  using  the  existing  procedure  for 
j P^(r)  P^(r)f(r)rdr  and  setting  P^(r)  “ 1,  f(r)  * 1*  Values  for  the 
integrals  of  0 can  be  obtained  with  the  existing  procedure  for  evaluating 
j DpQj(0)Qm(n)f (0)d0  by  setting  0^(0)  * 1,  f (©)  equal  to  the  expressions 
listed  above  for  f(0),  g(0),  and  h(0),  and  setting  Dp  - 1,  as 

needed.  Similarly,  values  for  the  z-integv:als  can  be  obtained  with  the 
existing  procedure  for  J Rk(z)Rn(z)dz  with  Rn(z)  • 1. 

The  effect  of  the  thermal  strains  is  to  change  the  right-hand 


side  of  the  system  of  equations  to  be  solved  by  subtracting  the  vector  of 
values  given  by  d/dF^"?  { JJJ  dV} . A subroutine  to  evaluate  the  thermal 

contribution  to  the  right-hand  side  of  the  equations  was  written  and  added 


to  the  program. 


APPENDIX  D 


THROUGH  CRACK  PROBLEM 


* mLm 


An  extension  to  the  analysis  was  made  to  handle  the  problem  of 
a through  crack  from  a cricular  hole.  Two  cracks  of  equal  length,  along 
a diameter  were  considered.  The  approach  taken  here  is  to  add  the  stress- free 
boundary  conditions  for  the  crack  faces  and  compute  an  energy  release  rate 
for  the  configuration.  The  requirements  for  stress- free  boundary  conditions 
are  given  in  Appendix  E.  The  energy  release  rate,  •* , is  the  amount  of  energy 
per  unit  of  length  that  the  laminated  plate  would  give  up  if  the  crack  were 
to  extend  a small  amount.  The  expression  for  is 


11m  ECC-^Q-E(C) 
AC*0  ^ 


(D-l) 


where  C is  the  crack  length,  AC  is  an  assumed  change  in  crack  length,  E(C+AC) 
and  E(C)  are  the  strain  energies  for  the  laminate  with  crack  lengths  C+AC 
and  C,  respectively,  and  t is  the  laminate  thickness. 

For  linear  elastic  materials,  the  strain  energy  is  equal  to  the 
complementary  energy  of  the  laminate.  The  expression  for  the  complementary 
energy  of  the  laminate  can  be  written  as 


>-v[cu  i 

p lc2i  c22 J *fp; 


(D-2) 


where  F is  the  vector  of  unknown  coefficients,  Fp  is  the  vector  of  prescribed 
coefficients,  and  [Cjj]  is  the  symmetric  flexibility  matrix.  Multiplying  out 
Equation  (D-2)  gives 


i r „t. 


^ I FiC,  ,F  + 2F*C12Fp  + F‘C22FpJ 


c 2 l nr 

The  solution  of  the  problem,  F*,  is  obtained  from  the  equation 


(D-3) 


Vc  * CllF*  + C12Fp  ‘ 0 


(D-4) 
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In  Equation  (D-4)  the  term  CI2Fp  18  the  right-hand  side,  RHS , of  the 
equations  to  be  solved.  Thus,  from  Equation  (D-4) 

CUF*  - -RHS  (D-5) 


Substituting  Equation  (D-5)  into  Equation  (D-3)  and  putting  the  solution 
F*  for  F gives 

"c  « | [ F*T(-RHS)  + 2F  T(RHS)  + FpC22Fp  j 
or 

nc  * 2 ! F*TrHS  + FPC22Fpj  (D-6) 


The  first  term  In  Equation  (D-6)  is  the  vector  product  of  the  solution  vector 
and  the  right-hand  side  vector.  The  second  term  In  Equation  (25)  Is 
contribution  due  to  the  prescribed  coefficients  Fp.  This  term  can  be 
evaluated  before  the  solution  vector  F*  Is  solved  for. 

Since  the  strain  energy  E is  equal  to  the  complementary  energy, 

^c,  Equation  (D-2)  can  be  written  as 


•b  * 11m 
AC"»0 


nc(C+AC)-nc(C) 

tAC 


This  procedure  for  evaluating  Jr  was  Included  in  the  program. 

A description  of  the  computer  program  and  the  input  data  is  given 
in  Appendix  F. 
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APPENDIX  E 


FORMULATION  OF  SEVEN  CRACK  PROPAGATION  PROBLEMS 
BY  IMPOSING  ADDITIONAL  STRESS  BOUNDARY  CONDITIONS 
ON  ELEMENTS  1,  3,  AND  5 


A crack  originating  at  the  inner  plate  radius  is  included  in  the 
finite  element  analysis  by  applying  the  boundary  conditions 


H-o-  ‘T«l 
e * iso0 


* T | * 0 

e*o°  0zle  - o° 

e * 180°  e » iso0 


to  elements  at  the  crack  location.  For  each  of  these  elements,  the  above 
boundary  conditions  reduce  to 


\ )i  Wik’A" + ruip;vJI 


9*0° 

0 - 180° 


(E-l) 


\ ] i {-Fuk<p;/r>'!JRk  + FCk(p*/t2,QJv|6 . 0.  • ° 


9 - 180° 


(E-2) 


SEE  ^{^(Pi/rJQ'R' 
i j k 3 3 K 


0° 

180' 


(E-3) 


Where 


p;  * 

pr- 


-P2 


-R2 


i //  _ r>  " 


-r; 


Q,l 

31  0 - O' 


9 - 180° 


q: 


Jl 


0° 

180” 


t | 1 odd  i values 
\ 0 even  j values 


f 1 even  j values 
l 0 odd  J values 
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(E-4) 


v V.r 


'4  m 


Equations  (E-l)  through  (E-4)  give  the  following  prescribed  F_k 
values  which  are  applied  to  each  element  at  a crack  location. 


F 

F 

F 

F 

F 


(1) 

ijk 

(1) 

ijl 

(3) 

ijk 

(3) 

ijk 

(3) 

Ijk 


F(0 

ij2  » 

0 . 


.(3) 

2ik 


k > 3,  all  i and  j 
all  i and  j 
all  j and  k,  i > 3 
i - 1,2,  j > 1 

all  k 
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APPENDIX  F 


INPUT  DATA  FOR  THREE-DIMENSIONAL 
EQUILIBRIUM  FINITE  ELEMENT  ANALYSIS 

The  analysis  consists  of  four  programs,  PROA,  PROB,  PROC,  and  PROD. 

PROA  read  in  the  data  to  describe  the  problem  and  generates  the  integrals 
of  the  products  of  polynomials  in  the  r and  z directions.  PRQA  also  integrates 
the  theta  modes  and  stores  them  on  tape.  There  is  an  option  not  to  evaluate 
these  integrals  if  they  are  already  evaluated  and  stored  on  the  tape.  PROA 
also  identifies  the  prescribed,  independent  and  dependent  coefficients  of 
the  stress  functions. 

PROB  has  no  card  input.  Data  is  transferred  to  PROB  from  the  tape. 

The  purpose  of  PROB  is  to  set  up  the  matrix  and  right-side  for  the  system 
of  equations  to  be  solved.  PROC  requires  data  to  describe  constraints  or 
independent  coefficients  that  are  associated  with  the  same  element.  These 
constraints  arise  from  satisfying  the  stress-free  boundary  conditions  at 
the  free  edge  around  the  circular  or  elliptical  hole.  Constraints  also 
arise  in  the  case  of  the  through-crack  problem.  The  form  of  the  constraints 
is  given  in  Appendices  B and  C.  The  purposes  of  PROC  is  to  solve  the  system 
of  equations  generated  by  PROB.  This  is  done  with  extended  core  storage 
or  with  mass  storage. 

PROD  evaluates  stresses  and  strains  at  selected  points  in  the  laminate. 
These  points  are  identified  by  input  data.  PROD  also  receives  input  data 
from  the  tape. 
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PROGRAM  PROA 


Variable 


JKEY,  NCP 


JKEY  * 1 for  temperature  problem. 


NCP  ' 8 n^crack0130^  j * temperature  problem  is  required,  Vead 

TEMPO  MlO.f 

TEMPO  * Temperature  change  \ 


TEMPO 

TEMPO  « Temperature  change 
and  ALP1 (3) , ALP2(3),  ALP(3) 

ALP1(3),  A2P2(3),  A2P(3)  are  the  coefficients  of  expansion 


1 for  circular  hole  problem.  Otherwise  a/b. 


3»3F10. 


Fia.9 


OPTQ  « 1 for  the  Q's  to  be  generated.  If  Q's  need  to  be  generated, 
read 


(3X,3F 


NTRE  is  the  size  of  the  Q matrix. 

RNX,  RNY,  RNXY . These  are  estimated  stresses, 


IPRNT,  ICHR,  ICHZ  312  ^ 

If  (IFRNT  .NE.O)  program  prints  check  information.  This  is 
used  for  checking  operation  of  program  and  transfer  of  data. 

ICHR  and  ICHZ  determines  whether  the  integrals  of  R and  Z 
should  be  recalculated.  Both  ICHR  and  ICHZ  should  have  the 
same  value.  If  a problem  with  the  same  element  R and  Z sizes  that 
are  on  Tape  1 is  being  run,  the  integrals  need  not  be  recalculated 
and  ICHR  ■ ICHZ  - 0.  Else  set  both  ■ 1. 

RB(l) , A (1)  2F10.0 

RB(2) , A(2)  2F10.0 

RB(l),  A(l)  - Inner  radius  and  R length  of  inner  elements 
respectively.  RB(2),  A(2)  - Irner  radius  and  R length  of 
outer  elements  respectively.  RB(2)  • RB(1)  + A(l) 

MMP(l),  MMP(2)  212 

ItiP(l)  and  1WP(2)  are  minus  the  starting  powers  of  R in  the  series 
for  the  inner  and  outer  elements  respectively.  Generally 
MMP(l)  » 2 and  MMP(2)  ■ 2 or  0 are  used  for  circular  holes. 


* i 


Var i able 


Format 


LNR  (set  equal  to  teio) 
NTZ  (set  equal  to  four) 


12 

12 

312 


MMR(l),  MMR(2) , MMR ( 3 ) 

MMR(l),  MMR(2)  and  MMR(3)  are  minus  the  starting  powers  of  the 
z-series  terms  going  from  the  center  to  the  outer  plies  respectively. 


LNRZ  (set  equal  to  zero) 

ZB(1),  11(1) 

ZB(2),  11(2) 


12 

2F10.0 

2F10.0 


ZB(3),  11(3)  2F10.0 

ZB  and  H are  the  starting  local  coordinates  for  the  z-axis  on  the 
element  and  the  thickness  of  each  ply  (starting  with  the  center 
ply ( 1 ) and  going  to  the  outer  ply(3))  respectively.  If 
MMR(I)  .CT.O,  then  ZB(I)  .GT.O. 

The  following  data  describes  the  theta  terms  in  the  series.  In  the 
following  description,  the  relation  between  N and  the  theta  terms  is 

N > 1 and  odd  N > 2 and  even 

Q(9)  ■ (cos  0)N  * Q(0)  « (cos  q)N  * sin  0 

These  data  are  read  in  for  each  stress  function,  and  for  the  theta 
variation  at  (1)  the  hole,  (2)  between  the  inner  and  outer  elements, 
and  (3)  at  the  outer  edge  of  the  outer  element. 

The  order  is 

Stress  Function  No.  I A . . . 

{0  variation 

0r  variation 
0rr  variation 

f 0 variation 

Location  B.  Between  inner  and  outer  elements  j 0r  variation 

1 0rr  variation 
. 0 variation 

Location  C.  Outer  edge  of  outer  element  j 0r  variation 

0rr  variation 


Stress  Function  No.  2 

Same  order  as  Stress  Function  No.  1 


Stress  Function  No.  3 
Same  order  as  above 
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Variable 


Forma  t 


For  each  stress  function  and  each  location  and  each  0,  0r  and  0rr,  two 
cards  are  required.  These  are 

NTV  12 

NTV  is  the  number  of  terms  in  the  theta  variation  series. 

NTV  < 18  and  (IDT  (ISF,I,J),  J - 1,  NTV).  2412 

IDT1  (ISF,I,J)  is  the  value  of  N described  in  the  preceding. 

ISF  * 1,  2,  or  3.  I goes  from  1 to  6 and  denotes  0,  0 , or  0 
variations.  r rr 

These  two  data  card  formats  will  be  repeated  a total  of  3ISF-S  x 
3 locations  x 3 (variations  on  0,  0r,  or  0rr)  - 27  times. 

The  next  set  of  data  describes  the  strain-stress  relations  and  the 
ply  singles.  For  each  ply,  starting  with  the  center  and  working  to 
the  outer  ply,  the  following  data  are  needed. 

PHI  F10.0 

PHI  is  the  ply  angle,  with  respect  to  the  9 equals  zero  line. 


PROGRAM  PROB  has  no  input  data. 
PROGRAM  PROC 


The  data  for  PROC  is  to  put  constraints  on  the  stress  function 
coefficients  that  could  not  easily  be  handled  elsewise.  These 
constraints  are  of  the  form,  * Yi  + Z^,  or  X^  * -y^.  Other 

constraints  of  the  form  Aj  * are  handled  automatically  within 

the  program.  The  data  read  in  are  (1)  the  number  of  equations  of 
the  form  X^  ■ + Z^  or  Xt  * -Yi,  and  (2)  identifiers  for  Xj, 

± *i,  «"d  ± zi*  The  input  is  as  follows. 

NEQ  12 

NEQ  is  the  number  of  equations  of  the  form  X£  ■ Y^  + Z^  or 
Xj  ■ to  be  satisfied. 

(J,DEP(K),  INI (K) , IN2(K),  K - I,  NEQ)  I2,8X,3(I6,4X) 

J is  not  used  by  program.  It  is  a means  to  number  the  input  data 

cards.  DF.P,  INI,  and  1N2  have  the  value  K + 10J  + 1001  + 1000ISF 

+ 10000  EL  where  the  element  number  is  EL,  the  stress  function 

number  is  ISF  and  the  coefficient  is  of  the  term  fijk  Pi(R)Qj(0)Rk(z)- 

A minus  sign  in  front  of  INI  or  IN2  indicates  there  is  a negative 
sign  in  the  equation.  DEP  is  always  positive.  INI  may  be  positive  or 
negative  but  not  zero.  IN2  may  be  positive,  negative  or  zero. 
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Variable 


Format 


r 


I 


PROGRAM  PROD 

IREAD,  I PUNCH  212 

IREAD  « 1 if  solution  vector  is  supplied  on  cards,  else 
IREAD  » 0.  IPUNCH  ■ 1 if  solution  is  to  be  punched  out  on 
cards,  else  IPUNCH  * 0. 

If  IREAD  - 1,  then  supply 

NC  110 

NC  ■«  Number  of  coefficients  in  solution 

(SV(I),  I - 1,  NC)  (2X.4E16.8) 

SV(I)  * solutions  vector. 

Comment:  The  following  is  data  describing  locations  in  the  laminate 


that  stresses  are  to  be  printed  out. 

L 12 

L ■ Element  No.  1 < L < 6.  If  L ■ 0,  program  stops. 

NR BAR  12 

NRBAR  m Number  of  radii  to  calculate  stresses  at  1 < NRBAR  < 11. 

(RBAR(I) , I - l,  NRBAR)  11F5.0 

NTBAR  12 

NTRAR  - Number  of  thetas  to  calculate  stresses  at  l < NTBAR  < 21. 

(TBAR(I),  I - 1,  NTBAR)  11F5.0 

Values  of  theta,  no  more  than  11  per  card. 

NZBAR  12 

NZBAR  * Number  of  Z locations  to  calculate  stresses  at  1 < NZBAR  < 11. 

(ZBAR(I),  I - 1,  NZBAR)  11F5.0 

Z locations. 


All  six  stresses  are  calculated  at  the  each  of  the  R,  9,  and  Z points 
described  by  input.  To  compute  stresses  for  another  set  of  points, 
repeat  input,  starting  with  L.  Last  card  in  data  for  PROD  should  be 
a blank  card. 


34 


REFERENCES 


(ll  Pipes,  R.  B.,  Kaminski,  B.  E.,  and  Pagano,  N.  J.,  "Influence  of  the 
Free  Edge  Upon  the  Strength  of  Angle-Ply  Laminates",  in  Analysis  of 
the  Test  Methods  for  High  Modulus  Fibers  and  Composites,  ASTM  STP  521, 

1973,  pp  218-228. 

(2]  Foye,  R.  L.,  and  Baker,  D.  J.,  "Design  of  Orthotropic  Laminates", 
presented  at  the  llth  Annual  AIAA  Structures,  Structural  Dynamics 
and  Materials  Conference,  Denver,  Colorado,  April  1970. 

(3l  Pipes,  R.  B.,  and  Pagano,  N.J.,  "Interlaminar  Stresses  in  Composites 

Under  Uniform  Axial  Extension",  Journal  of  Composite  Materials.  Vol.  4, 
October  1970,  pp  538-548. 

[4]  Rybicki,  E.  F.,  "Approximate  Three-Dimensional  Solutions  for  Symmetric 
Laminates  Under  Inplane  Loading",  Journal  of  Composite  Materials.  Vol.  5, 
July  1971,  pp  354-360. 

[5]  Pagano,  N.  J.,  "On  the  Calculation  of  Interlaminar  Normal  Stress  in 
Composite  Laminates",  Journal  of  Composite  “Materials.  Vol.  8,  January 

1974,  pp  89-105. 

[6]  Pagano,  N.  J.,  and  Pipes,  R.  B.,  "The  Influence  of  Stacking  Sequence 
on  Laminate  Strength",  Journal  of  Composite  Materials.  Vol.  5,  January 
1971,  pp  50-57. 

[7]  Pagano,  N.  J.,  and  Pipes,  R.  B.,  "Some  Observations  on  the  Interlaminar 
Strength  of  Composite  Materials",  International  Journal  of  Mechanical 
Sciences.  Vol.  15,  1973,  pp  679-688. 

[ 8 i Lackman,  L.  M.,  and  Pagano,  N.  J.,  "On  the  Prevention  of  Delamination 
in  Composite  Laminates",  presented  at  AIAA/ASME/SAE  15th  Structures, 
Structural  Dynamics  and  Materials  Conference,  Las  Vegas,  Nevada, 

April  17-19,  1974,  AIAA  Paper  No.  74-355. 

(9)  Dana,  J.  R.,  and  Barker,  R.  M.,  "Three-Dimensional  Analysis  for  the 
Stress  Distribution  Near  Circular  Holes  in  Laminated  Composites”, 

Report  VPI-E-74-18  from  Virginia  Polytechnic  Institute  and  State 
University  to  Department  of  Defense,  U.  S.  Army,  Contract  No.  DAA-F07- 
69-C-0444  with  Watervliet  arsenal,  Watervliet,  New  York,  August  1974. 

1 10 ] Rybicki,  E.  F.,  and  Hopper,  A.  T.,  "Analytical  Investigation  of  Stress 
Concentrations  Due  to  Holes  in  Fiber  Reinforced  Plastic  Laminated 
Plates,  Three-Dimensional  odeis",  Technical  Report  AFML-TR-73-100, 

June  1973. 

111]  Daniel,  I.  M.,  Rowlands,  R.  E.,  and  Whiteside,  J.  B.,  "Effects  of  Material 
and  Stacking  Sequence  on  Behavior  of  Composite  Plates  with  Holes", 
Experimental  Mechanics.  Vol.  14,  No.  I,  January  1974,  pp  1-9. 


35 


4 

[12]  Greszczuk,  L.  B.,  "Stress  Concentrations  and  Failure  Criteria  for 

Orthotropic  and  Anisotropic  Composite  Plates  with  Circular  Openings", 
Composite  Materials:  Testing  and  Design  (Second  Conference),  ASTM 

STP  497,  1972. 

[ 1 3 ] Waddoups,  M.  E.,  Eisenmann,  J.  R.,  and  Kaminski,  B.  E.,  "Macroscopic 
Fracture  Mechanics  of  Advanced  Composite  Materials",  Journal  of 
Composite  Materials.  Vol.  5,  October  1971,  p 446-454. 

[ 14]  Waszczuk,  J.  P.,  and  Cruse,  T.  A.,  "Failure  Mode  and  Strength  Predictions 
of  Anisotropic  Bolt  Bearing  Specimens",  AIAA  Paper  No.  71-354,  AIAA/ASME 
12th  Structural  Dynamics  and  Materials  Conference,  Anaheim,  California, 
April  19-21,  1971. 

[I5l  Rybicki,  E.  F.,  and  Hopper,  A.  T.,  "Transformations  for  Elliptical 

Holes  Applied  to  Finite  Element  Analysis",  Journal  of  Composite  Materials. 
Vol.  8,  January  1974. 


\ 


36 

& ‘OVfWWTNT  INC  III  : 


\ 


